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ABSTRACT 

We present spectrally resolved observations of triatomic carbon (C3) in several ro- vibrational transitions between the vibrational ground state and 
the low-energy v 2 bending mode at frequencies between 1654-1897 GHz along the sight-lines to the submillimeter continuum sources W31C and 
W49N, using Herschel's HIFI instrument. We detect C3 in absorption arising from the warm envelope surrounding the hot core, as indicated by 
the velocity peak position and shape of the line profile. The sensitivity does not allow to detect C3 absorption due to diffuse foreground clouds. 
From the column densities of the rotational levels in the vibrational ground state probed by the absorption we derive a rotation temperature (r rot ) 
of ~ 50-70 K, which is a good measure of the kinetic temperature of the absorbing gas, as radiative transitions within the vibrational ground state 
are forbidden. It is also in good agreement with the dust temperatures for W31C and W49N. Applying the partition function correction based 
on the derived r rot , we get column densities N(C 3 )~ 7-9xl0 14 cirr 2 and abundance x(C3)~10~ 8 with respect to H 2 . For W31C, using a radiative 
transfer model including far-infrared pumping by the dust continuum and a temperature gradient within the source along the line of sight we find 



that a model with x(C 3 )=10- 8 , r kin =30-50 K, N(C 3 )=1.5 10 1: 
agreement with the simple excitation analysis. 



cm fits the observations reasonably well and provides parameters in very good 
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1. Introduction 

Small carbon chains are relevant in the chemistry of stellar 
and interstellar environments for several reasons: ubiquitous in- 
terstellar spatial distribution, they likely participate in the for- 
mation of long carbon chain molecules, and they are prod- 
ucts in photo-fragmentation processes of larger species such as 
PAHs. Triatomic carbon, C3, was first te ntat ively identified in 
inters tellar gas bv lVan Orden et all dl995l) andHJaffner & Meverl 
(1995). The mid-infrared spectrum of C3 (V3 antisymmetric 
stretching mode) was measur ed in the cir c ums tellar envelope 
of CW Leo (IRC +10216) bv lHinkleet all (1 19881) . and in low- 
resolution interstellar ab sorption in the far- I R (V2 bending mode) 
toward Sgr B2 by ICernicharo et all d2000l) . lGiesen et all I (12001) 
discussed new laboratory data on the vibrational spectrum of 
C3 in its low-frequency bending mode and re-visited the first 
identification of the V2 R{2) line in absorption toward Sgr B2 
(IVan Orden et al.lll995h . The abundance and excitati on of C3 in 
transl ucent clo uds we r e dete rmine d convinci ngly by Maier et all 
d2001l) . lRoueffetal.1 d2002l) and lOka et alj d2003l) at optical 
wavelengths. 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



The Heterod y ne In strument for the Far-Infrared (HIFI; 
iDe Graauw et alj 20101) on bo ard the Herschel Space 
Observatory (Pilbr attet alJ [20 10), with its broad frequency 
coverage, high sensitivity and spectral resolution provides 
for the first time the opportunity for a systematic study of 
carbon chain molecules such as C3 through probing several 
ro-vibrational lines at full spectral resolution. In this Letter, 
we present the first results of our search for C3 obtained from 
observations of the sight-lines to the bright far-infrared (FIR) 
continuum sources W31C and W49N as part of the PRISMAS 
("PRobing InterStellar Molecules with Absorption line Studies) 
Key Program (Gerin et al. 2010). 

W3 1C (G10.6 - 0.4) is one of three bright HII regions within 
the W31 complex, and an extremely luminous s ubmillimeter 
and i nfrared continuum source (Lir ~ 1O 7 L Wright et al. 
119771) . Located at a distan ce of 4.8+°;^ kpc (Hj sh et al. 2003) 
with v LS r = -3 to -1 km s _1 dMiettinen et al J2 006). the sight-line 
to G10.6 - 0.4 intersects se veral foreground molecular clouds 
dCorbel & Eikenberry||2004l) . W49N is one of the three main 
IR peaks of W49A, which has a luminosity of (Lboi ~ 10 7 L 
Ward-Thompson & Robson 1990). It is located at a distance of 
11.4 kpc and has Vlsr = 12 km s ~ ' . The sight-line toward W49N 
is spectroscopically interesting because of the numerous features 
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Table 1. Spectroscopic parameters for the observed C3 transi- 
tions 



Name 


Transition 


Frequency " 


A-coeff 


E, 




(/',!/)<-« V) 


[MHz] 


10-v 1 


[K] 


P(10) 


(9,1)^(10,0) 


1654081.66(4.68) " 


2.38 


47.3 


P(4) 


(3,l)<-(4,0) 


1787890.57(6.90) 


2.72 


8.6 


2(2) 


(2,l)<-(2,0) 


1890558.06(0.25) 


7.51 


2.6 


2(4) 


(4,l)<-(4,0) 


1896706.56(0.15) 


7.58 


8.6 



R(0) 



" Experimental rest frequencies, uncertainties are given in parentheses. 
b Calculated frequency and \a uncertainty taken from CDMS catalog. 



has a perpendicular type spectr um with a calculat ed vibrational 
dipole moment of 0.437 Debye dJensen et al.l 19921) which shows 
prominent Q-, P-, and R -branch transitions. 

For the analysis of the HIFI C3 data, rest frequencies for 
^"(4), 2(2), and 2(4) were taken from laboratory measurements 
while for P(10) the rest frequency has been obtained from a 
global fit of all available laboratory data. 



Fig. 1. Energy level diagram of 2+ ground state and FI„ lowest 
bending mode state of C3. Due to nuclear spin statistics half of 
the rotational energy levels (dashed lines) are missing. Allowed 
ro-vibrational P-, Q-, and /^-branch transitions from v = to 
V - 1 follow - < — h selection rules. 



contributed by W49A itself, as well as by additional clouds as- 
sociated with the Sagittarius spiral arm (which crosses the line 
of sight twice). Numerous spectroscopic studies (both emission 
and absorption) hav e been carried out in the past toward s both 
these sources (e.g. JPlume et al]l2004UNeufeld et al.ll2010h . 



2. C , energy level diagram and radiative transitions 

Linear C3 has an energetically unusually low V2 bending 
mode at only ~ 90 K. Ro-vibrational transitions of the V2 
ba nd in its ^t ele ctronic grou nd state have been r eporte d 
by iGendrieschetalJ d2003l) and ISchmuttenmaer et all d 19901) . 
P(f) transitions with J = 2,4,6, R(J) with J = 0,2,4,6, 
and Q(J) with J - 2,4,. ..16 have been measured in the 
laboratory, using hi gh resolution Terahert z sideband spectrom- 
eters at Berkeley (Schmuttenmaeret al. 1990) and Cologne 



eters at Berkeley (^cnmuttenmaeret al. tyyu) ana cologne 
dGendriesch et al] 12003) with frequency accuracies of 7 MHz 
and 0.5 MHz respectively, which corresponds to line frequency 
uncertainties of 1.1 kms -1 and 0.08 kms _I . All data were 
used in a glo bal fit analysis to o btain most accurate molecu- 
lar constants dGiesen et alj 12001) which are presented as line 
lists for astronomical observations in the Cologne Data Base 
for Molecular Spectroscopy (CDMS) dMiiller et alJl2005l) . The 
ground state of C3 has £* symmetry while the vibrational ex- 
cited V2 is a two-fold degenerate bending state of n„ symmetry. 
Due to the 12 C nuclear spin of / = in the ground state only lev- 
els of (+) parity are present, while for the excited bending state 
only levels with (-) parity are allowed. As a consequence, in the 
ground state of C3 all odd numbered J rotational levels are miss- 
ing, whereas in the vibrational excited n„ state both, even and 
odd / rotational levels are present, but the two-fold degeneracy 
of the vibrational state is lifted (see Fig. Q]). Consequently, the 
statistical weights of the ro-vibrational levels is simply given by 
the rotational degeneracy: gj <v — 2 7+ 1. The V2 bending mode 



3. Observations and data reduction 

Along the sight-lines to W31C and W49N we have observed 
four lines of the v 2 bending mode, P(4), 2(2), 2(4) and P(10) 
(the latter only in W31C), of triatomic carbon in the upper side- 
band of the HIFI bands 7a, 7b and the lower sideband of band 
6b of the HIFI receiver. The observations of W31C and W49N 
were carried out on 2010 March 8 and 2010 April 19 respec- 
tively. The P(IQ) line was available as a "bonus" for an LO tun- 
ing dedicated to observe the CH line at 1661 GHz and it is yet 
to be observed in W49N with HIFI. All observations are in dual 
beam switch (DBS) mode and with the Wide Band Spectrometer 
with its spectral resolution of 1 . 1 MHz, corresponding to a ve- 
locity resolution of ~ 0.17 km s at the frequencies of the C3 
lines. To identify the line origin from the lower and upper side- 
bands, each line was observed with three LO settings shifted by 
15 km s _1 relative to each other. The 2(2) line also shows up in 
the 2(4) observations from the lower sideband, and, for one of 
the LO tunings, partially o verlap with the former. The data were 
first processed with HIPE dOtt et al .11201 Oh . and subsequently ex- 
ported to CLASS. At the high frequencies for these observations 
the H and V polarizations were at times found to be discrepant 
in the measured continuum level. Observations optimized for re- 
liable continuum measurement were used to select the spectra 
with the correct continuum level, used for the subsequent analy- 
sis. All spectra were smoothed to a resolution of ~ 0.68 kms -1 
and the rms noise level for the spectra lie between 0.01-0.03 K. 
Table [2] gives the measured double sideband continuum level 
(T c ). For the remainder of the paper we discuss the line inten- 
sities normalized to the single-sideband continuum level, where 
we have assumed a sideband gain ratio of unity. 

4. Results 

We have for the first time detected the spectrally resolved v 2 
band transitions P(4), P(10), 2(2) and 2(4) lines of C 3 . Fig. 
shows the observed spectra normalized to the (single sideband) 
continuum level. A multi-component Gaussian with common 
velocity width and spacing and individual amplitudes for each 
line was fitted to derive the basic parameters of the absorp- 
tion spectra. Table [2] presents the fit results and their uncertain- 
ties. The P(10) spectrum is affected by the spectral lines of CH 
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Fig. 2. The observed absorption spectra of the P{4), T^IO), 2(2) 
and 2(4) ro-vibrational transitions of C3 towards W49N {left) 
panel and W31C (right panel), plotted along with the best-fit 
profiles obtained by simultaneous Gaussian fitting to all lines. 



Table 2. Parameters derived from simultaneous Gaussian fit- 
ting to all the line profiles and column densities estimated 
from the fitted intensities. The fitted values of line center and 
linewidth are: for W31C V cen = -0.09+0.06 km s _1 and AV = 
4.5+O.lkms" 1 and for W49N V cen = 1 1.0+0.1 km s _1 , and AV = 
5.3±0.2kms _1 . 



( 1 66 1.1 0726 GHz) and H 2 0( 166 1.007637 GHz) from the upper 
sideband. 

The C3 absorption features towards W31C and W49N are 
centered near km s -1 and 1 1 km s _1 respectively. The systemic 
velocities of W31C and W49N are at Vlsr = — 1 kms -1 and 
12 kms -1 . Thus the C3 absorption lines detected here appear to 
be physically associated with the hot core itself and most likely 
arise in the lower density warm envelope surrounding them. We 
do not detect any absorption feature arising from the foreground 
clouds towards either of these sources (see Sec. [5}. 

4.1. Two-layer excitation analysis 

We first consider a simple two-layer model in which a (warm) 
absorbing layer without continuum opacity lies in front of a 
(hotter) emitting background source. We used the formalism ex- 
plained in detail in Appendix A to obtain estimates of the rota- 
tional temperature (r rot ) from the state specific column densities 
(Tab.|2). The J — 4 column densities are redundantly determined 
through both the P(4) and 2(4) transitions. One of the primary 
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Fig. 3. Model predictions for the Q(2), Q(4), P(2) and P(10) 
lines of the V2=l-0 transition of C3 in W31C. The black lines 
correspond to a model with n(H2)=10 5 cm 4 , x(C3)=5 10 s and 
7\in=50 K. The red lines correspond to a model with n(H2)=5 
10 5 cm 3 , ;t;(C3)=10 4 and 7^=50 K. Finally, the blue lines 
correspond to the same parameters as the red ones except that 
7W=30 K. 



sources of error in determination of the state specific column 
densities is the determination of the continuum levels. We esti- 
mate the uncertainty in the derived column densities to be 20%. 
For both W31C and W49N the discrepancy in the J - 4 column 
densities determined from P(4) and Q(4) is much larger than 
this uncertainty. We can only speculate on the reason: firstly, the 
assumptions of the simple two-layer model may not be appropri- 
ate and a more sophisticated model should be used (see below). 
Secondly, we note that in case the C3 absorption fills only a frac- 
tion of the solid angle of the continuum source within the beam, 
the intrinsic absorption would be much larger, possibly reaching 
saturation for the stronger 2(4) line. We use an average of the 
column densities at J - 4 obtained from P(4) and 2(4) for the 
rest of the paper. 

Using the formulae from Appendix A, we obtain for W49N a 
■Trot of 70 K with a large formal error. This is consisten t with 7dust 
derived from dust continuum observations as 43 K ( Vastel et al. 
2000) since r gas can be higher than 7d ust through direct gas heat- 
ing near the core. For T ml - 70 K we calculate from the 7=2 col- 
umn density a total column density of NXC3) = 9 10 14 cm" 2 . 

For W31C we have used the state-specific column densities 
(Nj) at 7=2, 4 and 10, to determine T rot and N(Ci) from the 



slope and the intercept respectively of a linear fit to the log( jjtt) 
vs Ej points ("rotation diagram"). Assuming a uniform 20% er- 
ror in the determination of the column densities we obtain T I0t = 
56+6 K and 7V(C 3 )= (7±0.5) 10 14 cnT 2 . The estimated T rot is 
consistent with the observed du st temperature of 52 K derived 
from continuum observations by Muell eret al.l (|2002). 

4.2. Radiative transfer models for C3 excitation in W31C 

In the likely scenario that the dust continuum source and its 
associated continuum opacity are coexistent with the gas that 
absorbs in the C3 lines, the C3 molecule will be embedded in 
a relatively strong continuum radiation which would contribute 
to the ro-vibrational excitation. In addition, the source intrinsic 
continuum will partially fill in the line absorption. Thus, in a 
more detailed approach we use a radiative transf er model, which 
considers FIR pumping by the dust continuum (Cerni charo et all 
2000) and a temperature gradient of the continuum source along 
the line of sight. For W31C we find that the C3 column densities 
can be interpreted by a cloud which is twice the size of the con- 
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tinuum source, has a molecular hydrogen density n(H2) = 10 5 
cm 4 , Vturb= 2 kms~', abundance x(C3) = 5 10~ 8 , kinetic tem- 
perature of 50 K and N(C 3 )= 1.5 10 15 cm 2 (see Fig. |3). The 
fact that the C3 column density comes out larger in this more de- 
tailed model is expected, as the dust opacity partially fills-in the 
absorption line. 

The main source of uncertainty in these models are the ro- 
vibrational collisional rates which are based on rather crude as- 
sumptions (ICernicharo et al.l l2000). However, irrespective of the 
adopted collisional rates the ground state is always thermalized, 
even in the presence of a strong IR radiation field in our models. 
We find that owing to the lack of a permanent electrical dipole 
moment the uncertainties in the collisional rates have little effect 
on the emerging intensities. Fig. [3] also shows that the resulting 
absorption depths are almost unaffected by a change in the local 
density (and hence abundance) by as much as a factor of 5, as 
long as the column density remains constant. We find that the 
V2=l-0 transitions are dominated by infrared pumping. The ro- 
vibrational excitation temperature in the inner part of the cloud 
is 35-37 K for the lowest rotational levels and around 20 for the 
higher-J states. In the external layers T ex decreases by 20 and 
10 K, respectively. In a second model with 11(^2) — 5 10 5 cm -3 
and jc(C3)=10~ 8 the results are almost equally consistent with the 
observations. We also note, that depending on the geometrical 
arrangement, the FIR-pumping in the ro-vibrational transitions 
can result in the net effect of lowering the rotational temperature 
in the vibrational ground state slightly below the kinetic temper- 
ature of the gas. 

We find that the major effect on the resulting absorption 
depths is related to the kinetic temperature adopted for the ab- 
sorbing gas. In Fig. [3] the blue lines correspond to a gas with 
7\in=30 K. The high-/ lines of the ground state are less popu- 
lated than in the previous case and the opacities for the V2 =1-0 
transitions decrease. If the opacity of the ro-vibrational lines be- 
come larger than 1, and the central continuum source is optically 
thick at the wavelengths of the V2=l-0 transitions, then the V2 
level becomes thermalized to the temperature of the dust. This 
effect can be counterbalanced by decreasing the gas temperature 
in the external layers of the cloud. Clearly, a more elaborate anal- 
ysis than can be presented in this letter, is needed to explore the 
full parameter range. 

5. Discussion 

In the absence of allowed radiative transitions in the ground 
state, the excitation of C3 in the V2 = state can be assumed to 
be thermalized to the kinetic temperature: T mt - 7\in- However, 
T I0t can be larger than r^st in the presence of direct gas heating 
mechanism like the photoelectric heating. The vi mode on the 
other hand could be excited by collisions and as well as by in- 
frared photons. With Einstein coefficients ranging between 2 and 
7 10 3 s , the line opacities can be high and the infrared pump- 
ing rather efficient. As a result, the excitation temperatures of the 
ro-vibrational lines are much lower, typically between the beam- 
diluted rd ust and 7\i n and hence the lines are seen in absorption. 
However, as explained above, the excitation temperature within 
the rotation ladder of the ground vibrational state can have val- 
ues larger than Tjust- 

Crude estimates based on the "two-layer" approach derive 
a T rot consistent with the dust temperature of the continuum 
source and give column densities of C3 in W31C and W49N to 
be between 7-9xl0 14 cirT 2 . Using the more elaborate radiative 
transfer model we derived a column density of 1.5 10 15 cm~ 2 for 
W31C for 7\i n of 30 K. Based on the discussion above we fur- 



ther argue that the thus derived column density of C3 is likely 
to be a lower limit. Implicit assumptions like a source filling 
factor of unity, r ex <sc T c etc. in the two-layer approach translate 
to an underestimate of T mt . Moreover the source intrinsic con- 
tinuum opacity at 1.9 THz partially re-fills the absorption (as 
shown in the radiative transfer model), and hence the column 
densities derived from the radiative transfer model are higher 
than those derived in the two-layer model. Thus taking all uncer- 
tainties into account we conclude that for both W3 1C and W49N 
the C3 column densities are ~ 10 15 cm 2 , correct to within a fac- 
tor of 2 or so. Th e H2 column densities in W3 1C and W49N are 
9.2 x 10 22 cirT 2 dMiettinen et al.ll2006l) and ~ 10 23 crrT 2 respec- 
tively, so that the abundances of C3 are ~ 10~ 8 . It is interesting to 
note that warm-up chemical models of the environment a round 
hot cores similar to the models by lHassel et al.l ([2008), with 
n-2 10 5 cirT 3 and Ay = 10 yield an abundance of 4-6xl0~ 8 . 

Based on absorption studies of C3 in optical wavelengths 
the C3 column densities in diffuse and tra nslucent clouds are 
found to range between 10 12 -10 13 crrT 2 dMaier et all 1200 U 
iRoueff et al.ll2002t lOka et al.ll2003l) . The C 3 column density ob- 
served in the present study is larger by about a factor of 100 or 
more than those from the optical studies. Thus, the non-detection 
of C3 in the foreground diffuse gas in the direction of our sources 
is consistent with the sensitivity of our observations. 
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Appendix A: Formulae used for the excitation 
analysis 

In the approximation of weak absorption (t « 1) and a negligi- 
ble population in the upper, V2 = 1, state, the lower state column 
density is given by: 



Ni = 



87TV 3 gl 



-f 



rdv 



The rotational temperature, r rot , 
specific column densities by 



(A.l) 



is calculated from the state 



TUU') 



In 



Nj./(2J'+1) 

N,/(2J+l) 



(A.2) 



where the energy of the levels is given by Ej - hBJ(J +1) 
and the rotational constant for the lower vibrational state ( vi = 
0)isB = 12908.242 MHz. 

Note that a ratio of level populations close to unity, 
( n iaj+\) ) ~ 1' i m pli es high values of T ml in comparison to the 
rotational energy scale defined by IiB/Icr, i.e. 0.62 K in the case 
of C3. The formal errors derived for the T vot are correspondingly 
very large. 

Assuming a thermalized population across the rotational lad- 
der with a unique value of T rot for all levels, a measured single 
state column density can be converted to the total column density 
N(Ci) using: 
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N, I hB 

N C , = P(T mt ) r-f- exp y—JiJ + 1) 

U + 1 \K-b1 rot 

T rot Nj ( hB 

exp\ ■ - J(J + 1) 



2hB/k B 27 + 1 



kBTrot 



(A3) 
(A.4) 



where the approximation T n 
the partition function 



» hB/k has been applied for 



P(T I0 t)= J] (2^+l)exp(-£ y /fer rot )^ 



7=0,2,4,. 



2hB ' 



(A.5) 



which is half of the usual value for a linear rotor due to the 
symmetry not allowing odd-J states. 
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